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ABSTRACT 

Context. Owing to their thermal emission, tidal disruption events (TDEs) were regularly detected in the soft X-rays and sometimes in 
the optical. Only a few TDEs have been detected at hard X-rays: two are high redshift beamed events, one of which occurred at the 
core of a nearby galaxy, and the most recent one is of a different nature, involving a compact object in the Milky Way. 

Aims. The aims of this work are to obtain a first sample of hard X-ray-selected unbeamed TDEs, to determine their frequency and to 
probe whether TDEs usually or exceptionally emit at hard X-ray energies. 

Methods. We performed extensive searches for hard X-ray flares af posifions in over 53 000 galaxies, up to a distance of 100 Mpc 
in the Swift BAT archive. Light curves were extracted and parametrized. The quiescent hard X-ray emission was used to exclude 
persistently active galactic nuclei. Significant flares from non-active galaxies were derived and checked for possible contamination. 
Resuits. We found a sample of nine TDE candidates, which translates into a rate of 2 x 10“^ galaxy^* yr“' above the BAT detection 
limit. This rate is consistent with those observed by XMM-Newton at soft X-rays and in the optical from SDSS observations, and is as 
expected from simulations. We conclude that hard X-ray emission should be ubiquitous in un-beamed TDEs and that electrons should 
be accelerated in their accretion flow. 
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1. Introduction 

Tidal disruption events (TDEs) occur when the gravitational po¬ 
tential of a compact object is large enough to tear apart an ap¬ 
proaching star or planet. The fate of the disrupted object depends 
on whether it approaches the tidal radius r, = (2 R*, 

where Mbh is the compact object mass and M, and R, are the 
mass and radius of the disrupting object. Part of the debris re¬ 
sulting from a disruption will fall on the compact object at a rate 
which depends on the disrupted object’s impact parameter and 
internal structure (binding energy). The accreted debris will heat 
up, might form a disc, accelerate particles in shocks, and possi¬ 
bly launch a jet. The heated material will shine in the optical and 
soft X-rays, accelerated particles may radiate in the hard X-rays 
through inverse Compton processes, and the jet, if aligned with 
the observer, may be detected as beamed emission and reach 
high apparent luminosities. 

The rate of TDEs at the centres of galaxies depends on the 
black hole mass and on the stellar population. In high-mass 
galaxies, N-body simulation (Brockamp et al. 20111 predicted 
rates of 3.0(+1.4) - 8.3(+4.2) x 10“^ galaxy ‘ yr ‘ for black hole 
masses in the range 10® - 10 ^Mq. 

Systematic search for bright X-ray flares in active and in¬ 
active galaxies using the Rosat All Sky Survey ( |Donley et al.| 
2002 I resulted in a rate of 1.8 x 10“® galaxy^^ yr"* for a peak 
flux detection limit < 10 erg s 'cm^ (unabsorbed). Excluding 
active galaxies, the rate dropped to 9.1 x 10“® galaxy^^ yr“'. Af¬ 
ter correcting for an estimate of the effect of local obscuration 
in the soft X-rays band, the corrected rate for TDEs increased to 
1.5 X 10“® galaxy^' yr“'. A slightly higher rate was obtained by 
[Khabibullin & Sazonov| ( |2014 1, who searched for fainter TDE 
candidates by correlating Rosat and XMM-Newton observations. 
The detection rate found in the XMM-Newton Slew Survey in 


non-a ctive galaxies is 2.3 x 10 ^ galaxy 
2008 1 , for a peak flux detection limit of 2.5 x 10 


-1 yr-I 


(Esquej et al. 


ergscm 


We note that signature for AGN activity was revealed by the [O 
III] 5007A/Hy6 line ratios in three out of the five TDE candidates 
detected in the XMM-Newton slew survey ( |Esquej et al.|2008} . 

Only two TDE candidates were detected in the optical from 
the SDSS Stripe 82 held (van Velzen et al. 201 l| l. van Velzen 


|& Earrar] ( |2014|l injected artificia l TDE flares (based on the 
simulations by Guillochon et al. 20141 in randomly selected 
SDSS targets and estimated an upper limit on the TDE rate of 
< 2 X 10“"^ galaxy^' yr“' (90% confidence level), in reasonable 
agreement with the rates determined from the X-ray observa¬ 
tions. 

The wide fields of view of the Burst Alert Telescope (BAT, 
2005|l on board the Swift satellite and of the 


Barthelmy et al. 


IBIS imager on board INTEGRAL ( [Ubertini et al.|[200^ allow 
a large fraction of the sky at hard X-rays to be monitored (the 
complete sky every few hours for 5w(/f/BAT). Two of the four 
TDEs detected at hard X-rays, Swift J164449.3-1-573451 ( Levan] 
let al.|201 11 IBurrows et al.|201 ![ jZauderer et al.|201 l|l and Swift 
J2058.4-H0516 ( |Cenko et al.|2012| l, are beamed, as confirmed by 
radio follow-up detections and luminosities that exceed the Ed¬ 
dington luminosity by large factors. The two other TDEs, discov¬ 
ered by INTEGRAL, were interpreted as the disruption of plane¬ 
tary objects, either by the super-massive black hole at the center 
of NGC 4845 ( Nikolajuk & Walter|2013 i or by a white dwarf in 
the globular cluster NGC 6388 ( Del Santo et al.|2014| l. 


2. Sample selection, data extraction and filtering 

Until now, non-relativistic candidate tidal disruption events with 
emission at hard X-rays were detected in nearby galaxies or in 
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globular clusters. As we aim to analyse SwiftfBAiT lightcurves for 
a large sample of candidate sources, we use the most complete 
catalogue of galaxies within 100 Mpc together with the list of 
globular clusters in the Milky Way. These objects were collected 


in the Gravitational Waves Galaxies Catalogue (GWGC, White 
et al.|20lT|l, which was built as an input for the search of grav 


itational waves from merging compact binaries and bursts from 
SNe, gamma-ray bursts, and magnetars. Included were 53 225 
galax ies with Lb > 10^ Lq and 150 globular clusters (Harris 
1996 1 . All TDEs that reach luminosities > 2.84 x 10"^"^ergT^ 


within 100 Mpc could, therefore, potentially be detected (as¬ 
suming a flare significance > 10 cr and duration > 10 days). 
This corresponds to the Eddington luminosity for a 2.3 x 10^ Mq 
black hole at 100 Mpc. 

The 5w//f/BAT reduction pipeline is described in |Tueller| 
et al.|(|20f0| and Baumgartner et al.](2013|l. Our pipeline is based 


on the BAT analysis software HEASOFT v 6.13. A first analysis 
was performed with the task batsurvey to create sky images 
in the eight standard energy bands using an input catalogue of 
86 bright sources (that have the potential to be detected in single 
pointings) for image cleaning. Background images were then de¬ 
rived, removing all detected excesses with the task batclean, 
and the weighted average determined on a daily basis. The vari¬ 
ability of the background was then smoothed pixel by pixel using 
a polynomial model with an order equal to the number of months 
in the data set. The BAT image analysis was then run again us¬ 
ing these averaged background maps. The image data were then 
stored in a database organised by sky pixel (using a fixed sky 
pixel grid with pixel < 5 arcmin) by accurately projecting the 
images on the sky pixels, whilst preserving fluxes. This database 
can then be used to build images, spectra, or lightcurves for any 
sky position. 

The result of our processing was compared to the standard 
results presented by the Swift team (lightcurves and spectra of 
bright sources from the SwiftfQPff 70-month survey catalogu^ 
and a very good agreement was found. The Swift team has cor¬ 
rected the spectral response derived from simulations to obtain 
a good fit of the very high signal-to-noise spectrum of the Crab 
nebula. The spectrum of the Crab nebula that we obtained shows 
similar deviations with respect to the simulations, but the devi¬ 
ations we derive are a factor of ~ 2 smaller than those obtained 
by the Swift team. We corrected the response in a similar way to 
the Swift team, but used ‘fudge factors’ that were a factor of ~ 2 
smaller. These modifications have no impact on the low signal- 
to-noise sources that we consider in the present study. 

5w//f/BAT images and lightcurves were built for every object 
of the catalogue from early 2005 to the end of 2013. Light-curves 
could be extracted for any time binning and energy range. For 
each time bin a weighted mosaic of the selected data is first pro¬ 
duced and the source flux is extracted, assuming a fixed source 
position and shape of the point spread function. The signal-to- 
noise of the source varies regularly because of its position in 
the BAT held of view with larger uncertainties observed every 
year. We tested various energy bands and Anally extracted the 
lightcurves that use the complete 20-195 keV band. This energy 
range was used throughout the whole analysis. Removing, for 
example, the highest energies did not increase the signal to noise 
significantly. 

As the X-ray emission of tidal disruption events is expected 
to last for several weeks, we decided to extract lightcurves with 
a resolution of five days, which is good enough to preserve the 


* http://swift.gsfc.nasa.gov/results/bs7Smon/ 
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general shape of the expected flares and allows a significant re¬ 
duction of the noise. 



Fig. 1. Distribution of positive and negative flares over significance- 
duration space in its most significant part. Red triangles show galaxies 
with positive flares, green stars shows star clusters and blue triangles 
show negative flares in the catalogue. 

Flare candidates were searched for using the following pro¬ 
cedure for each source. As a first step, we searched for the most 
significant bin in the lightcurve, then defined a flare duration by 
the period of time during which the flux remained above zero 
and calculated the global significance of the flare by integrating 
the lightcurve over that period. This was repeated for all excesses 
found in the lightcurve and an image of the most significant flare 
was finally built to measure its significance accurately. One flare 
was therefore kept for every source for further analysis. 

Because a number of the galaxies feature continuous X-ray 
activity, they are probably active galactic nuclei (or superim¬ 
posed with a galactic source). We excluded all sources from our 
flare candidates for which the time-averaged count rate, not in¬ 
cluding the most significant flare, was higher than 15% of the 
standard deviation of the count rate over the complete 5-day 
binned lightcurve. This criteria proved to select real persistent 
sources (or these contaminated by a nearby persistent source) 
effectively, with a total (integrated over eight years) significance 
> 4 cr (actually > 10 cr for most of them). This meant 288 
persistent sources were selected after the removal of confused 
flares (i.e. separated by less than 4 arcmin). Among them 115 
are known as AGNs, four as galaxies, seven as galaxy clusters 
and nine as galactic binaries, in the BAT catalogue ( |Baumgart^ 
|ner et al.|20131 l, where we use a maximum distance of 4 arcmin 
for source identification. In the remaining 153 persistent sources, 
we did not find a BAT source within 4 arcmin, however, for most 
of them a source of the BAT catalogue is closer to the fitted posi¬ 
tion than the original GWGC galaxy, indicating a clear contam¬ 
ination. In only 13 cases the persistent source had a PSF shifted 
by less than 6 arcmin from the GWGC galaxy and no counterpart 
in the BAT catalogue within 100 arcmin. These cases are further 
discussed in section 4.1. 

To understand the significance limit that should be applied 
to our list of flares, we used the same flare detection proce¬ 
dure as above on the negative lightcurves and compared the pos¬ 
itive and negative flare distributions in a significance-duration 
diagram. It turned out that flares with significance < 10 cr 
or duration < 40 days have a high probability to be spuri- 
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ous. Additionally we excluded flares within the ellipse matching 
((duration-40)/20)^-i-((significance-10)/5)^ < l)whereafew 
spurious negative flares could also be detected. Those regions 
of the significance-duration diagram were excluded to limit the 
chance of finding spurious flares in the positive flare sample. 

Finally we also rejected flares for which the integrated BAT 
images (about one degree in size) did show minima below -7 cr 
(or -1-7 cr for the negative flares) to eliminate high background 
events where count rates were not expected statistically. 

The distribution of the remaining positive and negative flares 
is presented in Fig.[T] Some negative flares below 15 cr and with 
a duration between 60 and 80 days are still detected, and should 
be taken into account in the analysis of the positive flares. 

We end up with 29 positive flare candidates (25 galaxies 
(two galaxy clusters) and four globular clusters) and eight neg¬ 
ative flares. The two galaxy clusters contained identified AGNs 
thus we removed those from our sample. Another two candidates 
correspond to the same flare, however their separation was high 
enough to avoid exclusion by the process described above. Thus 
we excluded the less significant of those two. At this point we 
are left with 26 candidates for further investigation. The positive 
flare parameters are listed in Table All fluxes and luminosities 
are given for the 20-195 keV energy band. The corresponding 
flare light curves and images are presented in Fig.[^ 

The point spread function of Swift BAT has a FWFIM of 17 
arcmin ( Barthelmy et al. [2005) 1. In some cases, several galaxies 
from our catalogue could therefore correspond to a single flare. 
In these cases, we selected as a counterpart the galaxy for which 
the flare appeared the most significant. These cases are further 
discussed in Section|3 

A visual inspection did not reveal any known X-ray emit¬ 
ting low redshift active galaxy in the vicinity of our flare candi¬ 
dates. In a few cases, discussed below, optically selected AGNs 
(Veron-Cetty & Veron 2010| l could potentially contaminate if 
they were to emit at hard X-rays. 


3. Flare candidate analysis 

We analysed each of the 26 flare candidates manually. As eight 
candidates were found in the ‘negative’ space, we expect some 
of the ‘positive’ candidates to be of instrumental origin. 

This verification started by checking the image of the field 
that was obtained during each flare. The point spread function 
was fitted to the image using a 2D Gaussian profile with a width 
cr = 8.2 arcmin (with mosaic_spec from the OSA toolkij^. 
The flare central coordinates are listed in columns PSF RA and 
PSF Dec of Table together with the positional uncertainty 
(Error column) at 90% confidence. The angular distance be¬ 
tween the most likely counterpart derived in Sect.j^and the cen¬ 
tre of the PSF is given in the column Offset. The error circles 
that result from the PSF fitting are indicated as the inner white 
dotted circles in Fig. while the outer white-black dashed cir¬ 
cles illustrate the radius of the PSF. When the most likely coun¬ 
terpart size is available it is indicated with a magenta circle in 
Fig. The difference between the counterparts and best fit po¬ 
sitions vary between 0.4 and 24 arcminutes. Differences larger 
than 10 arcminutes are most likely indicating a spurious asso¬ 
ciation. Very large (> 20 arcmin) differences indicate that the 
flare image is far from a point source and most likely related to 
instrumental or image cleaning problems. 

Also indicated in the flare images are the positions of other 
sources from the Gravitational Waves Galaxies Catalogue (red 


open circles) and from several other catalogues (the 13th edi¬ 
tion of the Catalogue of Quasars and Active Galactic Nuclei 
(cyan stars; [Veron-Cetty & Veron|2010i, the 7-year Swift-XRT 


point source catalogue (yellow squares; D’Elia et al.|2013|l, the 
Swift-BAiT survey of Galactic sources (green diamonds; Voss & 
Ajello|20T0)l, the Fe rmi LAT second source catalogue (blue tri¬ 
angles; Nolan et al.|2012|), the second Fermi/GBM GRB cata¬ 
logue (blue crosses; von Kienlin et al.||2014 , the list of super¬ 
novae observed by SwiftfKKY (cyan x; Russell & Immler|2012)l, 


the Swift X-ray observations of classical novae II (magenta x; 
Schwarz et al. 201 Ij l, the 5w(/f/BAT hard X-ray transient mon¬ 
itor catalogue (magenta stars; [Krimm et al. 201 3jl, the second 
Swift BAT GRB catalogue (red H-; Sakamoto et al.|20TT]l, the 70 
month 5wift -BAT all-sky hard X-ray survey (red x; |Baurngartner 
et al.|2013 1 obtained through the Vizier servic^ 


As very bright sources might affect the whole field of view of 
the BAT instrument and the source cleaning might be imperfect, 
we also checked whether our flare candidates occurred when any 
of the BAT bright sources were present within 60 degrees and 
flaring at the same time. This check was performed by compar¬ 
ing the flare candidate lightcurve with those of the very bright 
sources. Specific discussions are only required in a few cases 
and mentioned below. 

When a low significance flare occurs in regions with signifi¬ 
cant background, the fitting of the instrument point-spread func¬ 
tion to the flare excess can be unreliable and provide inaccu¬ 
rate results, in particular on the uncertainty of the flare position. 
These cases are discussed specifically. 

Swift observed gamma-ray burst afterglows for a long period 
of time with BAT pointing in a fixed direction and its field of 
view rotating on the sky because of Solar constraints. During 
this rotation, coding noise accumulated along arcs in the mosaic 
image at a level of typically (6 - 7) cr for an exposure of five 
days (see an example in Fig. [^. Such arcs are easy to spot, by 
analysing a wide field of view (at least three degrees), and some 
of our flare candidates do correspond to local maxima of an arc 
located by chance at the position of a GWGC object. These cases 
are discussed in the following subsections. 

In the following subsections, we describe each individual 
flare, sorted into four categories that depend on their most likely 
physical explanation. Only the sources listed in Sect. 3.4 are con¬ 
sidered as real TDE candidates. Eig. [^displays the field images 
and the BAT lightcurves, extracted at the positions of the flare, 
as derived from the point-spread function best fit. 


3.1. Galactic objects 

In this subsection we describe seven flare candidates, which 
likely originate in the Milky Way and were emitted by sources 
located either in globular clusters or in the Galactic plane. None 
of these sources are considered as a TDE candidate. 


3.1.1. Terzan 5 and IGR J17480-2446 

This flare is very significant (45.4 cr) and its image point-like. 
The flare onset occurred after 2010 Oct 01, peaked on 2010 Oct 
14, and probably lasted for about 72 days (up to 108 days). The 
flare, with a position consistent with the centre of the globu¬ 
lar cluster Terzan 5, was detected by INTEGRAL ( jBordas et al.j 
|2010| l and allowed the discovery of the new accreting pulsar IGR 
J17480-2446 ( Papitto et al.|201 f I. The second brightest flare in 
Terzan 5 peaked on 2012 July 18 and corresponds to a Type-I 


^ http://isdc.unige.ch/integral/analysis 


^ http://vizier.u-strasbg.fr/viz-bin/VizieR 
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Table 1. Flare properties and their host candidates. The flare’s properties (peak epoch, duration, peak flux, significance) were measured at the 
centre of the PSF fit, while the distance refers to the host itself and are taken from |White et al.H2011[ >. In the final column we list the most likely 
counterpart of the flare. The order of the table is as in Sect.|^ 


GWGC name 

Dl 

[Mpc] 

MJD 

Duration 

[days] 

Peak flux 
[10“'“ergs“^ cm“^] 

Significance 

Counterpart 

Galactic objects 

Terzan 5 

0.001 

55483 

72 

13.27 

45.4 

IGRJ17480-2446 

Liller 1 

0.001 

55812 

86 

2.46 

13.3 

MXB 1730-335 

NGC 6388 

0.001 

55771 

96 

3.23 

12.7 


NGC 6293 

0.001 

54583 

93 

3.52 

12.1 

RXJ1709.5-2639 

PGC 3079807 

65.2 

54798 

132 

25.37 

29.9 

Sw J1539.2-6227 

PGC 016029 

57.3 

55240 

95 

5.30 

17.0 

KTEri 

Gamma-ray bursts 

NGC 3835 

34.1 

53833 

47 

1.95 

17.1 

GRB 060319 

SDSS J114153.19H-510158.7 

97.5 

55125 

50 

3.68 

14.9 

GRB 091020A 

ESO 542-022 

76.0 

55471 

55 

2.91 

13.3 

GRB 100814A 

Relativistic TDE 

NGC 6213 

78.6 

55745 

121 

1.92 

15.2 

Sw J1644H-5734 

TDE candidates 

UGC 01791 

67.2 

53429 

82 

5.62 

13.7 


PGC 015259 

59.1 

55238 

60 

3.71 

10.4 


NGC 6021 

67.0 

53405 

71 

3.31 

10.1 


UGC 03317 

20.8 

55469 

61 

4.89 

9.8 


PGC 170392 

69.3 

53730 

69 

4.87 

10.6 


PGC 1185375 

22.9 

55258 

41 

6.69 

19.8 


PGC 133344 

97.1 

53480 

65 

3.08 

11.7 


PGC 1190358 

33.0 

55264 

108 

4.57 

10.4 


PGC 1127938 

76.8 

54116 

60 

3.57 

10.2 


Dubious cases 

NGC 1367 

17.3 

53737 

82 

3.99 

11.5 

SN 2005ke 

PGC 053636 

26.1 

53930 

81 

2.29 

10.4 


SDSS J114553.33H-635411.6 

58.0 

53828 

70 

1.54 

11.1 


NGC 6206 

72.5 

55750 

71 

3.22 

13.4 


ESO 161-001 

29.1 

54073 

90 

2.49 

11.3 


PGC 2464645 

31.2 

55846 

92 

2.70 

12.9 


PGC 3241255 

95.7 

55214 

65 

2.81 

9.9 



burst, which led to the discovery of the new accreting neutron 
star Terzan 5 X-3 ( |Bahramian et al.|2014| ). 


3.1.2. Lillet 1 


A significant flare was detected at the position of the globu¬ 
lar cluster Liller 1 from 2011 July 31 to 2011 Oct 25, peak¬ 
ing on 2011 Sept 07, with a total duration of approximately 86 
days. It was also detected by INTEGRAL and attributed to the 
rapid burster MXB 1730-335 ( Chenevez et al.|2011| l. The BAT 
flare image is centred 2.1 arcmin from the rapid burster. Flares 
from this source are regularly detected with 5w//f/BAT. The low 
mass X-ray binary 4U 1728-34, located 30 arcmin away, is much 
brighter but did not contaminate the flare significantly. 


3.1.3. NGC6388 

A significant flare was detected at the position of the globular 
cluster NGC 6388 from 2011 July 16 to 2011 Oct 20, peaking on 


2011 Sept 28, with a total duration of approximately 96 days. It 
was also detected by INTEGRAL and attributed to the disruption 
of an icy planet by a white dwarf ( [Del Santo et al.|[2014[ ). The 
low mass X-ray binary 4U 1735-44, located 31 arcmin away, is 
much brighter but did not contaminate the flare significantly. 


3.1.4. NGC 6293 and RX J1709.5-2639 

A significant flare was detected at the position of the globular 
cluster NGC 6293 from 2008 March 08 to 2008 June 09, peak¬ 
ing on 2008 April 26, with a total duration of approximately 93 
days. The BAT flare image is centred 8.8 arcmin away from the 
centre of NGC 6293, which is significant, but inside the globular 
cluster. The source shows recurrent flares in the BAT lightcurve, 
which are also detected by ISGRI and JEM-X on board INTE¬ 
GRAL. The most accurate source position, derived by accumulat¬ 
ing all JEM-X data, is RA = 257.378; DEC = -26.656; A = 0.7 
arcmin, which corresponds to that of the low-mass X-ray binary 
RX J1709.5-2639 ( |Degenaar etal.|2013^ 
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Table 2. Host name and position (Columns 1-5) and flare position (Columns 6-7) and uncertainty (Column 8) from the fit of the point spread 
function. The last column lists the differences between the GWGC candidate position and the PSF fit centre. 


GWGC Name 

Galactic 1 
[deg] 

Galactic b 
[deg] 

RA 

[deg] 

Dec 

[deg] 

PSERA 

[deg] 

PSE Dec 
[deg] 

Error 

[arcmin] 

Offset 

[arcmin] 

Galactic objects 

Terzan 5 

3.839 

1.687 

267.020 

-24.779 

267.038 

-24.782 

0.7 

1.0 

Tiller 1 

354.841 

-0.160 

263.352 

-33.389 

263.341 

-33.423 

1.9 

2.1 

NGC 6388 

345.556 

-6.737 

264.070 

-44.735 

264.095 

-44.691 

2.0 

2.8 

NGC 6293 

357.620 

7.834 

257.543 

-26.582 

257.398 

-26.650 

1.1 

8.8 

PGC 3079807 

321.085 

-5.726 

235.022 

-62.495 

234.797 

-62.470 

0.8 

6.4 

PGC 016029 

208.032 

-32.068 

71.949 

-10.235 

71.990 

-10.162 

1.3 

5.0 

Gamma-ray bursts 

NGC 3835 

137.527 

55.059 

176.020 

60.120 

175.946 

60.154 

1.3 

3.0 

SDSS J114153.19H-510158.7 

147.124 

62.703 

175.472 

51.033 

175.629 

51.071 

1.5 

6.4 

ESO 542-022 

169.255 

-76.835 

22.816 

-17.699 

22.615 

-17.880 

0.9 

15.8 

Relativistic TDE 

NGC 6213 

87.118 

39.815 

250.405 

57.815 

251.137 

57.732 

0.6 

23.9 

TDE candidates 

UGC 01791 

145.616 

-30.698 

34.973 

28.248 

34.977 

28.197 

1.7 

3.0 

PGC 015259 

199.738 

-33.581 

67.341 

-4.760 

67.344 

-4.816 

2.1 

3.4 

NGC 6021 

28.259 

45.576 

239.378 

15.956 

239.385 

15.957 

2.5 

0.4 

UGC 03317 

139.350 

20.771 

83.407 

73.724 

83.234 

73.688 

2.2 

3.6 

PGC 170392 

45.287 

-54.305 

336.693 

-15.023 

336.722 

-15.048 

2.7 

2.2 

PGC 1185375 

359.212 

48.971 

225.960 

1.127 

225.933 

1.091 

1.2 

2.7 

PGC 133344 

17.207 

-48.817 

325.733 

-30.133 

325.737 

-30.081 

2.4 

3.1 

PGC 1190358 

359.815 

48.778 

226.369 

1.292 

226.391 

1.303 

2.5 

1.5 

PGC 1127938 

138.276 

-63.103 

19.736 

-1.053 

19.681 

-1.137 

2.1 

6.0 

Dubious cases 

NGC 1367 

218.954 

-53.347 

53.755 

-24.933 

53.997 

-25.047 

1.6 

14.8 

PGC 053636 

358.020 

49.199 

225.263 

0.708 

225.416 

0.669 

1.8 

9.5 

SDSS J114553.33H-635411.6 

134.491 

51.729 

176.472 

63.903 

176.702 

63.832 

1.7 

7.4 

NGC 6206 

88.194 

39.831 

250.033 

58.617 

250.238 

58.585 

1.4 

6.7 

ESO 161-001 

266.355 

-26.674 

95.126 

-57.495 

94.956 

-57.548 

2.2 

6.3 

PGC 2464645 

101.430 

59.664 

211.342 

54.270 

211.232 

54.159 

1.7 

7.7 

PGC 3241255 

217.829 

-65.764 

40.067 

-26.550 

40.198 

-26.507 

2.2 

7.5 


3.1.5. PGC 3079807 and Swift J1539.2-6227 

A significant flare was detected at the position of the galaxy PGC 
3079807 from 2008 Nov 20 to 2009 April 01, peaking on 2008 
Nov 27, with a total duration of approximately 132 days. The 
BAT flare image is centred 6.4 arcmin away from PGC 3079807, 
which is significant and matches (within 0.02 arcmin) the posi¬ 
tion of LMXB Swift J1539.2-6227, marked with a magenta star 
in Fig. 1^ Our flare also matches in time with the discovery of 
Swift J1539.2-6227 ( |Krimm et all|20TT] l. It is not linked with 
PGC 3079807. 


3.1.6. PGC 016029 and KT Eri 

A significant flare was detected at the position of the galaxy PGC 
016029 from 2010 Jan 01 to 2010 April 06, peaking on 2010 Feb 
12, with a total duration of approximately 95 days. The BAT flare 
image is centred 5.0 arcmin away from PGC 016029, which is 
significant and matches (within 1.3 arcmin) the position of the 
classical novae KT Eri ([Schwarz et al.|[20TT]l, marked in Fig. 


with a magenta x. KT Eri erupted in November 2009 and became 
a super soft source in December ( [Bode et al.|2010 1 and reached 
an X-ray maximum in January 2010 ( Schwarz et al.|20rT] l. The 
hard X-ray flare detected by BAT is, therefore, not related to 
PGC 016029. 


3.2. Gamma-ray burst afterglows 

In this subsection, we list three flare candidates that correspond 
to known gamma-ray bursts. These gamma-ray bursts (GRB) are 
detected by chance and are not related to our sample of nearby 
galaxies. The BAT images of GRB are affected by circular arte¬ 
facts that result from the stacking of coding noise, which is gen¬ 
erated by many follow-up XRT observations that point to the 
GRB position, but rotated with respect to each other. In the cases 
reported here, XRT pointings took place for a significant amount 
of the flare duration and the BAT flare duration and fluxes are 
completely spurious. None of these sources are considered a 
TDE candidate. 
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Galactic Longitude 


Fig. 2. Illustration of the coding noise rings built during star observa¬ 
tions (here for an exposure of 121 days). This image is centred on GRB 
110328A, which is indicated by the red cross. 


3.2.1. NGC3835 and GRB 060319 

This flare, which occurred between 2006 March 16 and 2006 
April 30, featured a constant flux during the first 20 days and 
then faded away. Its total duration is between 47 and 135 days. 
The flare is shifted from NGC 3835 by 3.0 arcmin, which is out¬ 
side of the galaxy. 

The gamma-ray burst GRB 060319 ( [Sakamoto et al.|[2011| l 
occurred on 2006 March 19 for 8.9 s (shown as a vertical dashed 
line on the lightcurve in Fig. [^, when we detected a big rise in 
the BAT count rate. The GRB position is marked on the image 
in Fig. with a red cross which is within 15 arcmin of our flare 
position. 

3.2.2. SDSS J114153.19-1-510158.7 and GRB 091020 

This flare, which occurred between 2009 Sept 28 and 2009 Nov 
18, featured a maximum on 2009 Oct 10-21 (+2.5d) and decayed 
afterwards. Its total duration is between 20 and 50 days. The flare 
is shifted from SDSSJl 14153.19-H510158.7 by 6.4 arcmin. 

The gamma-ray burst GRB 091020A ( [Sakamoto et al.|2011[ 
[von Kienlin et al.|[2014| l occurred on 2009 Oct 20 (shown as a 
vertical dashed line on the lightcurve in Fig. |^, when we de¬ 
tected a rise of the BAT count rate. The GRB position is marked 
on the image in Fig. |^with a red cross, which is within 6.7 ar¬ 
cmin from our flare position. The difference between the flare 
and the GRB positions comes from the different periods of time 
used for the image integration. 


3.2.3. ESO 542-022 and GRB 100814A 

This flare, which occurred between 2010 Aug 10 and 2010 Oct 
04, has a rather undefined shape because of a limited signal to 
noise of 13 cr and a lack of features. Its total duration is ~ 55 
days. The flare is shifted from ESO 542-022 by 15.8 arcmin. 

The gamma-ray burst GRB 100814A (Gruber et al. 2011| l 
occurred on 2010 Aug 14 (shown as a vertical dashed line on 
the lightcurve in Fig. |^, which is consistent with our flare 

Article number, page 6 of[l4[ 


lightcurve. The GRB position is marked on the image in Fig. 
with a red cross, which is within 9.7 arcmin from our flare posi¬ 
tion. The difference between the flare and GRB positions comes 
from the different periods of time used for the image integration. 

3.3. Relativistic TDE candidates 

One of the relativistic TDE candidate detected by Swift occurred 
close to a GWGC galaxy and was therefore detected by our anal¬ 
ysis. 


3.3.1. NGC6213andGRB 110328A/Swift 
J164449.3-1-573451 


This flare, detected between 2011 March 26 and 2011 July 25, 
is very noisy and the image reveals that the emission is centred 
24 arcmin away from NGC 6213, on an unidentified persistent 
source (recognised as such by the criteria described in Sect. [^) 
The flare is not associated with NGC 6213. 

The gamma-ray burst GRB 110328A, which occurred on 
2011 March 28 (shown as a vertical dashed line on the lightcurve 
in Eig.[^, is consistent with the onset of the flare lightcurve. The 
GRB position determined by SwiftfXRT and indicated with a red 
cross in Eig.j^is within 9.3 + 0.6 arcmin from the flare’s best-fit 
position. The flare image shows extended background emission, 
related to the circular artefacts, and, therefore, the flare’s posi¬ 
tional uncertainty is underestimated and the lightcurve is uncer¬ 
tain. 

GRB 110328A/Swift J164449.3H-573451 is not an ordinary 


gamma-ray b urst. It was dubbed the longest GRB (Quataert & 
Kasen|20i2 i and was very bright (> 10"^^ ergs^'). XRT follow 


up revealed a new X-ray source but no UVOT detection was re¬ 
ported. This burst occurred in a small star-forming galaxy lo¬ 
cated at redshift z - 0.35 ( [Levan et al.[[201 1) . Its central black 
hole mass was constrained from X-ray and radio observations 
as log M = 5.5 + I.IMq Miller & Giiltekin] ( 201 l[ l. Several 


explanations were proposed to discuss the nature of this super 
long burst and, among them, a tidal disruption/fall-back accre¬ 
tion onto a massive black hole with collimated relativistic jet for- 
mation ([Burrows et al.|20lf Bloom et al.|2011[ [Zauderer et al.[ 


[201 1[ [Shao et al.|201 l[ ). [Krolik & Piran[P011| l addressed short¬ 
term variability during the X-ray flare by postulating a white 
dwarf disruption scenario, however the time-scale and fall-back 
rate analysis (Cannizzo et al. 201 1[ ) suggested a tidal obliteration 
event. Alternative mechanisms, such as the collapse of the galac¬ 
tic nucleus ( jPokuchae v & Erosh enko 2012| l or the emission from 
a millisecond pulsar ([Shao et al.|201 l|l, were also considered. 


3.4. TDE candidates 

This section presents flares for which we did not find any clear 
spurious counterpart or contaminant, and which are position- 
ally consistent with one of the GWGC galaxy. This sample 
brings together our best TDE candidates. TDE light curves are 
parametrized by a variability factor that is defined as the ratio 
between the luminosity of the flare’s peak Lpeak and the average 
luminosity over the quiescence period ^quiescence- 


3.4.1. UGC 01791 

This flare occurred between 2005 Jan 09 and 2005 April 01, with 
the most significant peak at 2005 Eeb 27. Its duration could be 
between 50 and 130 days. The flare is shifted from UGC 01791 
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by 3.0 + 1.7 arcmin and therefore, could come from that galaxy. 
We did not find any other source within the flare’s PSF radius. 
The closest XRT source is well outside of the PSF in the flare 
image. There is no circular artefact on a large-scale image. 

Assuming that UGC 01791 is the real counterpart, the peak 
luminosity and the energy emitted during the flare would be 3 x 
10"^“* erg s“* and 7 x 10^® ergs, respectively. The variability factor 
is > 47. 

3.4.2. PGC 015259 

This flare, detected between 2010 Feb 08 and 2010 April 09, 
features an initial increase with a slow decline over a period of 
60 days. The flare is shifted from PGC 015259 by 3.4 ±2.1 ar¬ 
cmin and could come from that galaxy. We did not And any other 
source within the flare PSF radius and there is no circular artefact 
on a large-scale image. 

Assuming that PGC 015259 is the real counterpart, the peak 
luminosity and the energy emitted during the flare would be 
1.55 X 10"^^ erg s“' and 3 x 10^° ergs, respectively. The variability 
factor is > 30. 

3.4.3. NGC 6021 

This flare was detected between 2005 Jan 02 and 2005 March 17 
and is featureless. Its duration could be 71 days but this is very 
approximate. The flare is shifted from NGC 6021 by 0.4 ± 2.5 
arcmin, and is, therefore, well centred on that galaxy. 

There is an other GWGC galaxy within the PSF radius, NGC 
6018, but we rejected it because it lies outside of the flare centre 
error circle and because the flare would be less significant for 
that position than our detection limit. There is no circular artefact 
on a large-scale image. 

Assuming that NGC 6021 is the real counterpart, the peak 
luminosity and the energy emitted during the flare would be 
1.78 X 10"^^ erg s“^ and 5 x 10^° ergs, respectively. The variability 
factor is > 28. 

3.4.4. UGC 03317 

This flare, detected between 2010 Sept 22 and 2010 Nov 22, 
features a maximum on 2010 Sept 29 and a decline that may last 
up to 61 days. The flare is shifted from the centre of UGC 03317 
by 3.6 + 2.3 arcmin. We did not And any other source within the 
flare’s PSF radius. There is no circular artefact on a large scale 
image. 

Assuming that UGC 03317 is the real counterpart, the peak 
luminosity and the energy emitted during the flare would be 
2.53 X 10"^^ erg s“^ and 4 x 10"^^ ergs, respectively. The variability 
factor is > 40. 

3.4.5. PGC 170392 

This flare, detected between 2005 Nov 28 and 2006 Jan 19, is 
quite noisy as the source reached the border of the BAT held of 
view during or after the flare. The flare position is shifted from 
PGC 170392 by 2.3 + 2.7 arcmin. We did not And any other 
source within the flare’s PSF radius and there is no circular arte¬ 
fact on a large-scale image. 

Assuming that PGC 170392 is the real counterpart, the peak 
luminosity and the energy emitted during the flare would be 2.8 x 
10"^“* erg s“' and 5 x 10^® ergs, respectively. The variability factor 
is > 35. 


3.4.6. PGC 1185375 

This flare, detected between 2010 Feb 03 and 2010 Feb 16, is 
highly significant and had a duration of 41 days. The significance 
maximum in the flare image is shifted from PGC 1185375’s cen¬ 
tre position by 2.7 ±1.2 arcmin. However, because the elongated 
background in this region significantly affects the result of the 
point-spread function fit, the GWGC galaxy and the flare could 
still be related. There is no clear contaminating source around, 
no point source detected by XRT, and no circular artefact on a 
large-scale image. 

There is an other GWGC galaxy at the border of the PSF 
radius, NGC 5831, but we rejected it as it lies even further out¬ 
side of the flare centre error circle and because the flare would 
be less significant for that position. AGNs lying within the PSF 
radius are too distant (z > 0.85) to contribute to this emission 
effectively. The elongated structure is probably instrumental. 

If we assume that the flare was emitted by PGC 1185375, the 
peak luminosity and the energy emitted during the flare would be 
4.2 X 10"^^ erg s“* and 10^® ergs. The variability factor is > 51.5. 

3.4.7. PGC 133344 

This flare, detected between 2005 March 17 and 2005 May 21, 
lasted for about 65 days. The significance maximum in the flare 
image is shifted from PGC 133344’s centre position by 3.1 ± 2.4 
arcmin. There is a possible second close-by weak BAT source in 
the held (at more negative latitude and similar longitude, about 
20 arcmin from the flare position) that could easily explain this 
shift in position. 

If we assume that the flare was emitted by PGC 133344, the 
peak luminosity and the energy emitted during the flare would 
be 3.5 X 10"^^ergs“^ and 1 x 10^* ergs. The variability factor is 
> 14. 

The second GWGC galaxy within the PSF radius is less sig¬ 
nificant (9 vs 11 sigma), more distant and further away from the 
PSF centre, and shows a shorter flare. AGNs within the PSF cir¬ 
cle are at high redshifts (z>1.6). We note that the optically iden¬ 
tified quasars 2QZ J214250-2946, at redshift 2.8, lies 1.2 arcmin 
away from our flare candidate and PGC 133344 (|Veron-Cetty &| 
|Veron|20101 l. 

3.4.8. PGC 1190358 

This flare, detected between 2009 Dec 29 and 2010 April 16, 
lasted for about 108 days with a clear maximum on 2010 March 
08. 

The significance in the flare image matches PGC 1190358 
with a distance of 1.5 ± 2.5 arcmin. The extended image does 
not show circular artefacts. 

If we assume that the flare was emitted by PGC 1190358, the 
peak luminosity and the energy emitted during the flare would be 
6 X 10"^^ ergs^^ and 2 x 10^“ ergs. The variability factor is > 34. 

3.4.9. PGC 1127938 

This flare, detected between 2007 Jan 08 and 2007 March 09, 
lasted for about 60 days with a peak on 2007 Jan 15. 

The significance maximum in the flare image is shifted from 
PGC 1127938’s centre position by 6.0 ±2.1 arcmin. There is no 
contaminant nor circular artefact in the extended image. 

If we assume that the flare was emitted by PGC 1127938, the 
peak luminosity and the energy emitted during the flare would be 
2.5 X 10"^"^ erg s“' and 4x 10^*’ ergs. The variability factor is > 29. 

Article number, page 7 of 


14 







A&A proofs: manuscript no. tdeflares 


3.5. Dubious cases 


3.5.5. ESO 161-001 


Here we list the remaining flares. They have a low signal to 
noise, and the flare images do not show a point-like excess. None 
of these sources are considered to be TDE candidates. 


3.5.1. NGC 1367 


This flare, detected between 2006 Sept 17 and 2006 Dec 16, 
lasted for about 90 days. 

The flares does not look at all like a point source and is an 
artefact, which results from XRT pointings on GRB 060729, lo¬ 
cated 4.8 degrees away. 


A significant flare was detected close to the galaxy NGC 1367 
from 2005 Nov 08 to 2006 Jan 29, peaking on 2006 Jan 01, with 
a total duration of approximately 82 days. The BAT flare image 
is centred 14.8 arcmin away from NGC 1367, so does not match 
with the galaxy. 

The gamma-ray burst GRB 090909B ( |von Kienlin et ah 
2014| l occurred in the same region as our BAT flare but occurred 
four years later and so does not match. 


The supernovae la SN 2005ke which exploded on 2005 Nov 
06 + 2 days (with optical and X-ray maxima occurring ~ 20 and 
respectively ~ 45 days later pmmler et al.||2006 i) is 14 arcmin 
away from the BAT flare but these two events coincide in time. 
The large scale image is affected by a circular artefact centred on 
SN 2005ke, we must conclude that this flare is an artefact from 
the observations of SN 2005ke. 


3.5.6. PGC 2464645 

This flare, detected between 2011 Aug 19 and 2011 Nov 19, 
lasted for about 92 days and featured a peak on 2011 Oct 
12. It was found to match with PGC 2464645 and SDSS 
J 140357.26+542228.6 (these were not associated with the dedu¬ 
plication algorithm, which uses a smaller distance as criteria). 

The significance maximum in the image is shifted from 
PGC 2464645’s centre position by 7.7 ±1.7 arcmin. Near to 
PGC 2464645, many XRT pointing observations was conducted 
within MlOl galaxy probably as a result of SN 201 Ife follow¬ 
up. SN 201 Ife was discovered on 24 Aug 2011 in MlOl at a 
position shifted by 17.5 arcminutes from the detected flare posi¬ 
tion. 

3.5.7. PGC 3241255 


3.5.2. PGC 053636 

This flare, detected between 2006 July 11 and 2006 Sept 30, 
features an early maximum and a slow decrease over 81 days. 
The lightcurves shows very large error bars after the flare when 
the source moved close to the border of the BAT field of view 
and a gap near the beginning because two months of data were 
missing from our analysis for technical reasons. 

The flare is shifted from PGC 053636 by 9.5 ±1.8 arcmin 
so they are not related to each other. There are several optically 
selected AGNs close by but no known X-ray source in the XRT 
catalogue. The flare around PGC 053636 originates from a large 
ring-like structure (like the one shown in Fig. |^, which is cen¬ 
tred on GRB 060814, and which is completely unrelated and 
separated by more than 20 degrees from PGC 053636. 


3.5.3. SDSS J114553.33+635411.6 

This flare, detected between 2006 Feb 19 and 2006 April 30, fea¬ 
tures a rather flat lightcurve during about 70 days and a rather 
low significance. The flare could be longer as it is immediately 
followed by a period of two months during which data were 
missing in our analysis for technical reasons. The flare is shifted 
from SDSS J114553.33+635411.6 by 7.4 ± 1.7 arcmin. The ex¬ 
tended image shows circular structures, the flare is an artefact re¬ 
sulting from the follow-up observations of GRB 060319 located 
3.9 degrees away. 


3.5.4. NGC 6206 

This flare, detected between 2011 May 11 and 2011 July 20, 
lasted for about 71 days and featured a peak on 2011 July 07. 

The significance maximum in the image is shifted from NGC 
6206’s centre position by 6.7 ±1.4 arcmin. Circular artefacts 
indicate that this flare resulted from XRT follow-up observations 
of GRB 110328A/SwiftJ164449.3+573451, located 1.2 degrees 
away. 


This flare, detected between 2009 Nov 21 and 2010 Jan 25, 
lasted for about 65 days. The significance image is elongated 
and does not correspond to a point source. This flare is an arte¬ 
fact that results from follow-up observations of GRB 091127, 
located 8.3 degrees away. 


4. Discussion 


In this work we focused on the search for flares produced by 
TDEs that originate in non-active galaxies. It was, thus, impor¬ 
tant to identify and remove persistent sources. As a result we 
found five new persistent sources, which were not known as ac¬ 
tive galaxies previously and discuss them in Sect. 4.1 The four 
flares that originate from globular clusters are discussed in Sect. 


Most flares were detected close to non-active GWGC galax¬ 
ies. Among them, three correspond to gamma-ray burst after¬ 
glows that are located by chance close to NGC 3835, SDSS 
J114153.19+510158.7, and FSO 542-022. The nine remaining 
TDF candidates are discussed in Sect. 14.31 


4.1. New persistent hard X-ray sources 

The removal of persistent sources is described in Sec. where 
persistence is defined as when the averaged count rate (exclud¬ 
ing the most significant flare) divided by its standard deviation 
is larger than 0.15. In this section we take a closer look at the 
153 persistent sources not previously known as active galactic 
nuclei. For the majority of them we found that a source of the 
BAT catalogue was closer to the X-ray source PSF centre than 
the GWGC galaxy. In only 13 cases the persistent source had a 
PSF shifted by less than six arcmin and no counterpart in the 
BAT catalogue within 100 arcmin. The reason that we found 
new sources not listed in |B aumgartner et ^ ( |20 1 3] l could be ex¬ 
plained by the longer time period for the data in our analysis, 
meaning higher significance for persistent sources. We checked 
other X-ray catalogues but did not find any obvious contaminant. 
In addition, we checked for contaminants in the 13th edition of 
the Quasars and Active Galactic Nuclei Catalogue (|Veron-Cetty| 
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|& Veron||20101l. In five cases, our targets are close enough to 
known optical AGNs (UGC 3255, Zw 022.021, ESO 500-G34, 
ESO 439-G09, and UM 163), which are not listed in the BAT 
catalogue. 

We considered the remaining eight persistent sources as 
AGN candidates and rejected four of them as being dubious or 
contaminated for the following reasons. Several circular arte¬ 
facts, as presented in Fig.[^ affect the lightcurve of UGC 512 and 
mean that this AGN is probably a spurious candidate. As PGC 
2793378 was considered a candidate for an obscured galaxy 
cluster ( |Woudt et al.|1999) l and as a Fermi AGN, 2EGL J1629.6- 
6141, lies 4 arcmin away, this AGN candidate is also dubi¬ 
ous. Finally the galaxies SDSS J21 1238.66 -hI 12822.8 and PGC 
022813 should also be disregarded as they may be contami¬ 
nated by the AGN UGC 11700 and by the galaxy cluster RX 
J0808.7-1-6731, respectively. 


We investigated the light curves of the four remaining AGN 
candidates using different time bins and coding fractions. Two 
sources (ESO 385-014 and PGC 070775) could be excluded as 
their light curves were affected by coding noise and showed en¬ 
hancements only at specific time periods. Two persistent objects 
remained: ESO 315-020 and ESO 153-003. The significance of 
the corresponding excesses in the total BAT image is 13.2 and 
4.8, respectively. Thus the last one is more uncertain, however 
an XRT point source was also detected near the position of ESO 
315-020. We conclude that these two GWGC galaxies are likely 
to be new active galactic nuclei. 


The study of ESO 315-020 is difficult because it is contami¬ 
nated by a bright foreground star ( |de Grijs & Peletier|2000 1 and 
because it is considered as an edge-on spiral galaxy (|Kregel et al.| 


2002 20051, which makes the detection of activity more diffi¬ 


cult. Optical observations of ESO 153-003 revealed a dust lane 
along its major axis which, together with a radial velocity mea¬ 
surement along the minor axis, allowed for its classification as 
an edge-on SO galaxy (Hau & Forbes|2006 1. 


4.2. Globular cluster 

We also detected four flares at the positions of globular clus¬ 
ters included in the GWGC catalogue. The most significant flare 
was detected in Terzan 5 and most probably originated from the 
LMXB/NS EXO 1745-248. The second flare, detected in Filler 
1, could be related to the rapid burster MXB 1730-335 ( |Homer| 
|et al.|20()T] ). The third source is the TDE candidate that was dis¬ 
covered by |Del Santo et aL] ( |2014| l and modelled as the disrup¬ 
tion of an icy planet by a white dwarf. Our last candidate orig¬ 
inated within the angular size of the cluster NGC 6293 and is 
most probably associated with the low mass X-ray binary RX 
J1709.5-2639. If this association is correct, the flare luminosity 
would be 4.2 x 10^^ erg s“'. 


4.3. TDE candidates 

Our sample contained 26 flare candidates close to GWGC galax¬ 
ies. After excluding all cases showing quite clear indications of 
contamination, such as presence of GRB afterglows, supernova, 
classical nova. Galactic binaries, and dubious sources, we are 
left with nine candidates. These flares occurred near the galax¬ 
ies UGC 01791, PGC 015259, NGC 6021, UGC 03317, PGC 
170392, PGC 1185375, PGC 133344, PGC 1190358, and PGC 
1127938. 



Lp^ak [ srg S 


Fig. 3. Histogram of flare’s peak luminosity in the 5-days binned BAT 
light curves in the TDE candidates sample. 


If we assume that all these nine flares originated from tidal 
disruption events, the TDE rate, above the BAT detection limit, 
could be estimated as 


Ntdf — 


9 events 


52874 galaxies x 8 years 


2.1 X 10 ^ yr 'galaxy 


( 1 ) 


This should formally be considered as an upper limit, as 
some of these flares could have a different origin. 

The peak luminosities of the flares were estimated using 
the peak fluxes derived from the 5-day binned BAT count-rate 
lightcurves, converted to physical fluxes using the Crab as a 
template, as well as the distances to the host galaxies. The his¬ 
togram of these luminosities, presented in Fig. spans values 
JQ43-44 gj.g expected for black holes of 10^^ Mq. 

The duration of the flares observed by BAT are between 40 
and 110 days and their energy is in the range 10^"“^* ergs. These 
are of course lower limits since TDEs are expected to last longer, 
below the BAT sensitivity. These durations are approximate and 
have a tendency to increase with the maximum flux of the flares, 
which qualitatively match the expectations. 

During the six-month exposure of the Rosat All-Sky Survey, 
IDonley et al.| (|2002|| detected two TDEs from galaxies of the 
GWGC (NGC 5905 and IC 3599). Our hard X-ray TDE rate 
is, therefore, comparable to what was derived at soft X-rays. A 
more accurate comparison, taking into account the instrumental 
sensitivity limits and the GWGC sample is required. To derive 
the TDE distribution expected in a sample like the GWGC, we 
need to estimate the black hole masses in the host galaxies, the 
mass function of the disrupted objects, and the distribution of 
impact parameters. 

A black-hole mass was estimated for each galaxy of the 
GWGC. We decided to use the total galaxy luminosity as a num¬ 
ber of them lack morphological information and to account for 
irregular galaxies. The absolute B-band luminosities were con¬ 
verted to black hole masses following [McConnell & Ma] ( |2013| l. 
For the distribution of the disrupted objects masses M*, we used 
the Kroupa initial mass function (as parametrized in Equation 


20 of|Stone & Metzger|2014||. The impact parameter distribu¬ 


tions were derived for each black hole and disrupted object mass 
assuming isotropic encounters. 
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Fig. 4. Expected distribution of peak fluxes for TDE in the GWGC for 
all black holes (blue) and only from those with masses in the range 10®“’ 
Mq (green). The red histogram shows the TDE candidates detected by 
us. The vertical line shows the expected eRosita flare sensitivity at hard 
X-rays. 


We defined a grid of parameter space with 5<log Mbh< 10 
(step of 0.1), -l.l<log Mt <1 (step of 0.01), and inverse of im¬ 
pact parameter -0.2<log (J3 = rtlrp)<Q.6 (step of 0.01). The grid 
was limited in the region where the impact parameter is larger 
than the black hole gravitational radius. On every point of that 
grid, we computed the peak accretion rate of a TDE, using equa¬ 
tion A1 from |Guillochon & Ramirez-Ruiz| ( [2013| l, and converted 
it to a peak bolometric luminosity, using the radiative efficiency 
given by Fig. 1 in |Ulmer| ( |1999| l, and finally to a hard X-ray peak 
luminosity assuming that 4% of the total emission is released at 
hard X-rays, as in active galactic nuclei ( Risaliti & Elvis|2004) l. 

Each point of the grid was weighted by the Kroupa IMF and 
by the probability of the impact parameter. The peak luminosity 
was then summed in the grid, using the black hole mass distri¬ 
bution derived from the GWGC as a function of the distance to 
finally obtain the peak flux distribution. This distribution was 
renormalised such that the number of events predicted above 
2 X 10“'® erg s“'cm“^ (our BAT sensitivity limit) corresponds 
to the observed hard X-ray TDE rate. 

The resulting distribution is shown in the blue histogram 
in Fig. 1^ which shows two components corresponding to the 
change of equation of state between low and high mass stars. 
The green histogram shows the contribution of black hole masses 
of 10®“’ Mq. The red histogram is derived from the observed 
sample. The shape of the red histogram does not reproduce the 
high flux tail of the blue histogram but we should remember that 
above a flux of 10“® erg s“'cm“’, the blue histogram predicts 
one TDE every five years for the GWGC, which is consistent 
with our non-detection. The blue distribution is also based on 
several assumptions (like black-hole mass and IMF) that may 
not be representative for the GWGC sample, for instance three 
of our nine TDE candidates occurred in irregular galaxies. 

Some interesting numbers can be derived from Fig. 1^ and 
be compared with previous TDE rate estimates. The TDE rate, 
which is derived from the blue histogram above a flux of 
2.5 X 10“'^ ergs“* cm“’ is 1.3 x 10“"' ga laxy ' vr“'. matches well 


the rate of 2.3 x 10 "'galaxy ' yr ' ([Esquej et al.||2008|l de¬ 


rived from the XMM-Newton slew survey (and optical flares as 
well). The TDE rate, for black holes in the range 10®“’ Mq, 
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Fig. 5. Average spectrum of the nine BAT TDEs. 


which was obtained by integrating the green histogram, is 4.9 x 
10“^ galaxy“'yr“', which also matches the expectations from 
the simulations of |Brockamp et al.| ( |201 l) l well. The TDE rate de¬ 
rived from the blue histogram above a flux of 10“'’ erg s ' cm“’ 
(the Rosat flux limit) is 1.2 x 10“"'galaxy“' yr“', which is sig¬ 
nificantly higher than the observed rate of 3 x 10“® galaxy ' yr“' 
( |Khabibullin & Sazonov|[2()14[ ). As the latter is also in tension 
with the XMM-Newton observations, it is difficult to conclude if 
this has a specific meaning. 

Starting in 2016, eRosita ( [Merloni et al.|2012l l will scan the 
X-ray sky regularly, observing any sky position every six months 
or less. Taking into account that one TDE flare could be detected 
for three months, we expect seven TDE per year to be detected at 
hard X-rays for GWGC galaxies and more from higher redshift 
sources. Most TDEs that originate in the GWGC will be detected 
at soft X-rays, which will give an observational access to most 
of the histogram featured in Fig. 

We have built the average spectrum of our nine TDE can¬ 
didates by extracting the spectrum of each individual flare and 
averaging them. The average spectrum, unfolded using a pow- 
erlaw model, is shown in Fig. The average spectrum is well 
fitted by a power law with an index of 2.18 + 0.08 and an aver¬ 
age flux of7.8xl0“" ergs“' cm“’ (20-100 keV). The power-law 
slope is similar to that observed in active galactic nuclei. 


5. Conclusions 

We have detected nine TDE candidates at hard X-rays during 
eight years in a sample of ~ 53 000 galaxies closer to 100 
Mpc. These events are not beamed as their peak luminosities 
(3 - 34 X 10"'’ ergs“') are not significantly higher than the Ed¬ 
dington limits (10"'’ ^“"'® erg s“') derived from the corresponding 
black hole masses. 

Assuming that 4% of the TDE bolometric luminosity is emit¬ 
ted at hard X-rays, the observed hard X-ray fluxes and rate are 
consistent with the rates observed by XMM-Newton at soft X- 
rays or in the optical from SDSS observations. It is also consis¬ 
tent with expectations from simulations. 

We conclude that hard X-ray emission should be ubiquitous 
in unbeamed TDEs. The average hard X-ray spectral shape of 
our TDE candidates is, in addition, very similar to that usually 
observed in active galactic nuclei, which points towards inverse 
Compton processes and electrons that are accelerated in shocks 
in the accretion flow. We also discovered a long flare of RX 
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J1709.5-2639 and two new persistent active galactic nuclei can¬ 
didates. 
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Fig. 6. Left panels: the light curves derived at the position of the flares maximum, with red individual points with error bars, while the detected 
flare is shaded with red. If a GRB was detected nearby, a target galaxy vertical dashed line shows the time of the GRB. Right panels: significance 
images built over the flare period in logarithmic greyscale (the lower limit is fixed at 2 cr). The associated GWGC galaxies are marked with 
red-filled circles, while other galaxies from the catalogue are shown with empty red circles. Cyan stars mark AGNs, yellow squares are XRT 
point sources, green diamonds show BAT galactic sources, blue triangles are Fermi sources, blue and red + show GRBs from Fermi and BAT, 
cyan x shows supemovae candidates, magenta x mark classical novae detections, red x code bright BAT sources, and magenta stars present BAT 
transients. The big black empty circle represents PSF fit with its radius equal sigma of the Gaussian (8.2 arcsec) with inner circle with radius equal 
to the error of the PSF. Magenta circles shows galaxy/globular cluster size if it was given in the catalogue. Overplotted contours are significance 
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Fig. § (continued) 





Article number, page 13 of 


14 













































































































































































A&A proofs: manuscript no. tdeflares 



apnjijen 3U3eie0 




lO 




in 

o 

U 

y T 


- ' __ 

:— 


' - 

' ^ j 


a 

i 






Fig. § (continued) 
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